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ABSTRACT 



Design of linear control systems is examined using 
algebraic methods. Previous work indicates that the system 
describing equations are nonlinear. Various combinations of 
feedback and cascade compensated systems are analyzed 
using velocity constant, bandwidth and complex root speci- 
fications, Methods for !i n'eari z i ng^ these system equations 
are demonstrated and examples are solved using specific 
digital computer programs. Feasibility of a general computer 
program is discussed. 
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I . I ntroduct i on 

Design is a word of many meanings and forms. Generally, 
it means a method or procedure for determining the value of 
system parameters such that the system wi I I respond in a 
prescribed manner. These methods are usually derived for 
synthesis in the time domain by applying modern state space 
techniques or in the complex frequency domain using root 
locus. Bode, Nyquist, or Nichol’s plot. 

All the methods have their advantages and disadvantages. 
Time domain analysis with state space allows the designer 
to examine the entire system at one time. He can only select 
one set of parameter values and investigate their effect on 
the system. The advent of digital computers greatly aided 
this process because the computer is easily adapted to a 
systematic trial and error design procedure. 

The complex frequency plane methods are inherently 
trial and error. These methods are graphical and extremely 
time consuming in their application. Present day programming 
is gradually coming to the rescue but it has limitations in 
the amount of useful and legible information that can be 
displayed on a graph. 

An outgrowth of the graphical methods was Mftrovic’s 
method which was generalized by the Parameter Plane technique. 
This procedure transforms the analytical model of the system 
into a graphical representation. The benefit of this method 
is that it can handle mu 1 1 i parameters . 

Q 



From the parameter plane technique system design has 
returned to the analytical world and developed a multi- 
parameter procedure called algebraic methods. Algebraic 
methods reduces the system differential equations and spec- 
ifications prescribing the desired response into a set of 
simultaneous algebraic equations. The simplicity of the 
method is that it finds the parameter values that collect- 
ively provide the required system response. The method has 
developed over the years but its implementation on computers 
has not been explored. It is the intention of this thesis 
to investigate this area. 
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1 1 . Bas i c Theory 



The backbone of the algebraic methods is expressing the 

systems characteristic equation as a polynomial of the form 

TL < 

Z a*S - O (2-1) 

where s is the comp I ex var i ab I e and is defined as 

$- - S UJn 4- | — 5* ' (2-2) 

and the a^'s are real coefficients which are functions of 
the system parameters. Two independent equations in a k ,S , 
and w n are obtained if equation (2-2) is substituted in 
equation (2-1) and the real and imaginary parts set equal 
to zero. Performing additional manipulation and introducing 



a new variable (1), (2), we obtain 

i 0 .^ 10 * &_,($)= o 

K~o 




(2-3) 


d K u)ia 0(t.(S) — O 

V • A 




(2-4) 


lx — U 

where 


2(S)J 


(2-5) 


with 

0.(S)-o, 0,{S) - -1 amo 


VUS)=-&($) 


(2-6) 



Other values ofp^jfor various values of S are tabulated 
in Appendix I. Observe that for a real root S= - cr and 
substituting in equation (2-1) gives 



wr vr 

E (-rf«u*=o 

K-o 



(2-7) 



Previous work provided graphical methods for portraying 
root locations as a function of any two system parameters. 

(2), (3), (4), (5). If we are willing to abandon graphical 
analysis then the equations can include any number of unknown 
parameters. Proceed algebraically by substituting root 
locations in the basic equations (2-3) and (2-4) giving two 
equations as functions of a I I the unknowns. Requirement for 
solution is the number of equations equal the number of un- 
knowns. One otL-ii.s idditional equations by one of two 
philosophies; first, substitution of additional root locations 
into the basic equations; second, derivation of additional 
algebraic relationships from system performance specifications. 
Determination and solution of the proper number of equations 
results in a simultaneous set of parameter values that satisfy 
stated requirements. Naturally, if the parameter values are 
not satisfactory the performance specifications must be altered 
or the additional equations generated changed and a new 
solution affected . 

Inherently these equations provide a set of nonlinear 
simultaneous equations to be solved. The non I i neari ti es are 

generated from the system configuration or conversion of the 
performance specifications into algebraic relationships. The 
purpose of this thesis is to develop methods of solving these 
nonlinear equations. Basic approach will be that of linear- 
ization by substitution or manipulating the form of system 
components to facilitate linearization. When linearization 
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fails it is hoped the system will lend itself to one of the 
standard analytical methods for solving linear simultaneous 
equat ions (6 ) . 



III. Third Order Development 

The basic system general block diagram is shown in 
Figure 3-1. In this section we will use a type one, third 
order plant to demonstrate the different solution techniques. 
Naturally, it is required that some type of compensation be 
employed, namely, feedback, cascade and feedback-cascade 
combi nati on, so that the system satisfies the specifications. 
The third order system allows for development of the method 
while keeping the equation length reasonable. As needed 
reference will be made to higher order systems and some 
higher order examples will be solved. 

3.1 Feedback Compensation 

Consideration will first be given to tachometer feed- 
back, Figure 3-2, where the system gain (K) and the tachometer 
gain ( K -^ ) are unknown. Solution for the unknowns requires two 
simultaneous equations. Since the most common system speci- 
fication is a pair of dominant complex roots, equation (2-3) 
and (2-4) will be implemented. 

The system characteristic equation is 

S i + tPi+Pi)S iL + (PlP2 + KKt)S+ MO (3-1) 

from which the a^ coefficients and the given S and are 
substituted in the root equations giving 

P2+ PZ)W^($)P ( 3-2 ) 

*.<k(S)i-(piP2+ 4>,lSH (H+ P2)0)^(S)fW,?<4(S) = 0 (3-3) 
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It can be seen these equations are nonlinear due to the 



KK-t term. Linear equations are obtained by making the 
following transformation of variables 



After making the substitution and rearranging we get the 
following matrix form 







X 




'-f 1P2W„A(S)-(PUV2>)^15)- U)^ <MS) ' 


6M 




X 




-P1P2 -W l <t>S\ 



Linear equations are easily solved by digital comouters 
with programs using matrix algebra for implementing the 
standard linear equation solving techniques. It is a simple 
operation to inverse transform the linearizing substitution 
and obtain the parameter values 



Observe that for this specific case the matrix equation 
degenerates to give 



X- K 

Y- K&t 



(3-4) 



vc~x 



and 




(3-6) 






9 -, is) 



(3-7) 



and 



v- - Pifie #,cs H j> 3 (s) 

because 0 O (5) = 0 



(3-8) 
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Another possible method of feedback is acceleration 
feedback. Figure 3-3, which is very similar to tachometer 
feedback but affects the second order coefficient of the 
characteristic equation instead of the first order coeffic- 
ient. This type of compensation generates a system with two 
unknowns, therefore, requiring two equations which we will 
obtain from the root specifications. The system character- 
istic equat i on i s 

sVcPl+fi+KJOsS P1PZS+K. - o (3-9) 



Appropriate substitution in the root equations yields 



K£,(S> Wi DVWS)+ (pHPE+-KK4Ia)*0,19 +<*)* fe(S) -O (3-io) 



and 



R9fc(S) 1-PI P2 VOr,^,(S) + ( P1+ PZ+ - c> 



(3-11 ) 



As in the tachometer feedback case the following linearizing 
transformation of variables can be made 



Making the substitution and rearranging gives the matrix 
simultaneous equations 



( 3 - 12 ) 



'0,(5) W»0,($) ‘ 


"X" 


\ 


-PLP2UXv<KC5)-(.Pi+F2)t) * <#>CS) -W * <4(S) 


_&(5) wi&<s) 


x 




•pi ?i w^fpivK.) 



(3-13) 
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For this type of feedback we do not get the degeneracy 
due to the r ortunate location of 0 O (S)» but it is seen the 
matrix equations for tachometer and acceleration feedback 
are the same except for the second column of the coefficient 
matrix. Therefore, we can expect that when both types of 
feedback are used simultaneously the matrix equation will 
be very similiar to (3-5) and (3-12). At this point we will 
state the results if a fourth order system is used. The 
matrix equation is 



w'fysj 

wo wi«u& 



KKa 



J L 



“«Pi?2Wn.4^'(PVP2 + ?l?3^P3)W^4(S)'(Pltp^B)w^<()»iS) 

Wrv A &) '(P1P2 tPiBf PiPi)U)!v4‘ 

(3-14) 



From this it is seen that the form is the same as the third 
order case and all the increase in system order did was add 
a term to the right hand side. 

Let us now combine the types of feedback, Figure 3-4, 
and solve the problem. Immediately it is seen that before 
we can proceed we have to develop another equation besides 
the root equations. A frequently used specification is 
error constants which has different meaning depending on 
the type of system. For the type one system this means we 
can use the velocity constant defined as 



Kv - SGoCs) 



(3-15) 
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with the steady state error for a ramo inDut, R(s) 

Css = - U M _L_ a _L 

s - ?0 u-^CS) *->o SQo(s) 

The characteristic equation for the system is 

sV(Pi-fl?Z + K^)sV (PlPZ-tKW) StKiO 



I /s 2 



( 3 - 16 ) 

(3-17) 



and the forward transfer function is 

6 .(s)= t 

5 [ t‘i+?l)fe+?2) +-KKi\| 

Ado lying the root equations and the velocity coefficient 



(3-18) 



equation we obtain the following set of simultaneous non- 



linear equations in K, KK-|., and KK a 

Kv?lPZ + KvKK*-K=.0 (3 _ |g) 

0-, i?)K+ (pipz-v <kis) 4> iCi) -a (3 _ 20) 

<^C?)K+C? 1 P 2 + (3 _ 2 , } 



As in the seoarate type of feedback compensation we can apply 
the same transformation of variables to linearize the equations 

Xc-K 

V- \06fc (3-22) 

2 = 



Making the substitution and rearranging we obtain the matrix 
equation 



1 

o 

> 

V 

U 


X " 




--Kvf>L?2 


<k?) W^(S)|»^K) 


r 


- 


- P1P2 W*. 4> 0 (s) -C Pl+P2)0) 1 <j> t t?) - W *<£»(?) 




* 

L , 




_ 
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As predicted the above result is very similar to the 
feedback cases examined separately. The coefficient matrix 
second column is from the tachometer equation (3-5) and the 
acceleration equation (3-13). Similarly, the right side is 
not changed at all. From this we conclude the type of 
compensation for the root location specification affects terms 
in the coefficient matrix while the right hand side is affected 
only by the system order. 

A possible specification for generating an extra equation 
is bandwidth. The conventional definition of bandwidth will 
be used, that is, the frequency at which the magnitude of the 
output is one half the input magnitude for the closed loop 
transmission with S = jW^. Because the magnitude involves 
squaring the transfer function we can expect to gain a 
nonlinear equation of a rather difficult form. For feedback 
compensation, it is found the problem is not extremely difficult. 

Consider the tachometer feedback. Figure 3-2, where the 
open loop transfer function is 

— - (3-24) 

and the closed loop transfer function is 

QcCS)= — (3-25) 

S(Sv?iXsirVZ)-+- KKtS 

to which we apply the bandwidth requirement giving 

1 - - h ( 3 - 26 ) 

2. - [f>lp2-uo + [K- W b z (m P2}] 

1 9 



To simplify the exoression make the substitution 



A- -U)b a (Pl+?Z) 



and 



B = Pi? 2. - tOb 2- 



(3-27) 



whereby clearing gives 



- K? + 2A K - 2E .UV- VCKt ■+ A*- ft* - uO b 2 K*Kt - O 



(3-28) 



Before proceding any further let us examine the second equation 
required for solution. Select the velocity constant specifi- 
cation because a complex root specification generates two 
equations and it is not our intention to deal with over 
determined systems, hence using 



Examination of equations (3-28) and (3-29) shows an extreme 
nonlinearity. Fortunately, if we consider KK-j- as an unknown 
we can use the elimination method to reduce the system to a 
quadratic equation in forward gain K. Solving equation (3-29) 
for KK t 



Kv - L'wi sf. 
s-p-o L 




gives 



-VCvPlVi 



(3-29) 



K Kt = ( X- KvW?2)/Kv 



(3-30) 



and substituting this value in equation (3-28) gives 



[.-UO-UJb/Kv 1 ] VC + L2A - ZU) J/Kv ZWbPiPZ/W} K 



-vlWfc?! Pz-v C PlPzf ] - o 



(3-31 ) 
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This is easily solved for the system gain K yielding two 
values which can be substituted in equation ( 3 - 30 ) to give 



two values of «-)- 




(3-32) 



At this point the designer has a choice of which set of gains 
to use for the solution. Usually it is not a difficult choice 
because one of the gain values is negative. Situations where 
both gains are positive values the designer can pick the values 
that suit his purpose. 

Now we’ll shift our attention to the acceleration feed- 
back, Figure 3-3. From previous expressions it can be expected 
the feedback gain in the second order term will help the 
solution. As in the tachometer case we will use bandwidth and 
velocity constant specifications. The velocity coefficient is 



Kv - K-/P1.P2. 



giving one equation 



K.- Kv'PIVZ. 



(3-33) 



The bandwidth specification gives 



1 - 
2 . “ 



- + ZT-fttfK&u + £ - Bj z 



(3-34) 



where 



A= P1P2-6J6 7 - 



B- LUi?(pl + Pl) 






Rear rangement and substitution for K from (3-33) gives a 



a quadrat i c in K a . 
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(PI + (-1 tOt ( PIP 2. PZ K*) P*- 

+ te x -coZ~/o' 1 - zbPi pz Pv - (Pi PZ kv) Z J -- o 



(3-35) 

Again this gives two values for the acceleration and forward 
gains. As in the tachometer case the choice is left to the 
designer but usually one of the choices will have an unde- 
sirable negative gain. 

3.2 Cascade Compensation 

Cascade compensation is a widely used method of compen- 
sation. In general, it is commonly referred to as a lead, 
lag, or lead-lag filter. In this development we will not pre- 
specify the type of filter or a criteria for realizability. 

The decision of realizability is left to the system or filter 
d e s i g n e r . 

The first case examined is the single section compen- 
sator, Figure 3-5. In keeping with the major requirement we 
need three equations to solve for the unknown system gain (K), 
filter zero(Z), and the filter pole(P). If given a complex 
root specification and a velocity constant we can generate 
the required equations. A possible unknown is eliminated by 
considering any gain associated with the filter to be included 
in the overall system gain. 

The system characteristic equation is 



S*-h(Pl+PZ-hP)$ + ( ?(P±+Pl)i-PJ.PZ ) sV (&P2 P-+K)5~tK£~0 



(3-36) 
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indicating one nonlinearity, the KZ term. The velocity 
constant specification yields 



k 



V - 



KLs+i) 

(s +?)(3l-rd)Cs+ Pi) 



S -O 



KZ 
PIP IP 



( 3 - 37 ) 



which has the same nonlinear term as the characteristic 
equation, making the now familiar transformation of variables 



X - Kb 



( 3 - 38 ) 



provides linearity with the three unknowns P, K, X which after 
affecting the root equations gives the matrix equation 



o -KS&J.P2. 




' K 


1 


o 


tOnfa PJ.P2 (P1 +PZ)UJa 




p 


~ j 


! -FlPZtiu.fa -(pj+n?) 










u)^4i - tv,? 


UJ't 4, Pl?2Li),\ l <t>o 




K 


1 

1 

! 


1 

- PJ PI Ultfa -IPUK2.) 








i 


fa - fa ( 3- 



After solution we once again inverse transform the variable 
f or the filter pole 



£= X/K 



( 3 - 40 ) 



Unlike the feedback system a change in the system order 
affects not only the right side of the matrix equation but, 
also, the left hand coefficient matrix. The reason for this 
effect can be seen if the reduction of the system to a forward 
transfer function with unity feedback is examined. The filter 
pole multiplies all poles of the basic plant, thus inserting 
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the filter pole in all the characteristic coefficients 
except a n and a 0 . 

The next logical consideration is a multisection cascaded 
filter. To demonstrate the best analysis form we will use the 
double section filter, Figure 3-6. The addition of each filter 
section adds two more unknowns, the section zero and pole, 
which requires two extra equations per section added. At the 
risk of completely specifying the system roots the easiest 
specification to use giving two equations is a root specifi- 
cation, otherwise other specifications will have to be trans- 
lated into algebraic equations. The root specification can 
be very helpful since if not given explicitly it can be selected 
such that it does not effect dominance of any specified roots. 

Before proceding further, to ease the notation, multiply 
the plant denominator giving 

f , in 



where 

A= P.'+P-J 
6= P.'Vx 



(3-42) 



Reducing the overall system by normal methods, we achieve the 

following character i st i c equation 

(A+Pl + P2)s^ + ( &+&(Pl + PZ) + P±P2.)s -+(B(PJ+Pi) 

+*(PJLPZ) *K)$ 

At this point, note the transformation of the unknowns into 
the seven groupings PI, P2, PIP2, KZ I , KZ2, K, and KZIZ2. This 
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double section may not be that unwieldy to solve but the 
situation gets worse as the sections increase. A simple 
manipulation can alleviate this undesirable condition. 

Conversion of the filter to polynomial form vice factored 
form prevents the generation of undesirable nonlinear terms. 
Granted some remain but they can be handled by appropriate 
variable transformation. The polynomial form of the filter 
i s 



X$ V Y 
+ 2 



(3-44) 



where X - 2.i -t ~2z. 

Y- 

P,-h Pz. 

?,Pz. 



(3-45) 



Using this filter form in the system reduction we achieve a 
new and simpler characteristic equation 

s 6 -t tliA/i- 2 ) 7 *- CBVJ V A2 f k) S 

(32“ XX ) 5 &Y — O 



The nonlinear terms can be linearized by 

v5= XX 



(3-47) 



The application of the root equations twice and the velocity 
constant 



Kv « Lim 5 fc(sV X5~f Y) ^ _KY 

s ~*° (s x -hWs-t2)(sh/i &-£ 



(3-48) 
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yields the matrix simultaneous equation form 



K 



[R 1 




Da 



where for root specification one the coefficients 



AO, O'- <&(*) 

A (i, 2 .) ~ Aw* 4,(0 4 Z ( s) y Bat* 4>o(t) 



A(h 4 )~ 

A( h s) ~ / A) CU,t 4>z( ~A 4.z ( %) 

A ( z ,0 - &*($) 

A ( Z , l ) “ 4 , ($) ($) ( s ) 

AU t l ) ~ 

A&A)- 

A (*,*)- &uJ^4> t ($) +tu H 4 <fa(s) 

SCO ~ -B0Jn4i(S) -Arif fo(f) - M if Ms) 
8 (z) *■ - B w<4 (S) ~A lA Aft) - c*>£t 4>s($) 



are functions of zeta one and omega one. For root specifi- 



cation two the third and fourth rows have the same form as 



rows one and two, ie, 

/jfaj c)— A(2,L ) 



and 



&(l) ” BO) 
BU) ~ &iz) 



- 49 ) 



-50) 



- 51 ) 
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where A(3,i), A(4,i), B ( 3 ) and B(4) are functions of zeta two 

and omega two. The remaining row is 

A(Sj) “ -/• O 

A(s,’l) -- AV £ 

(3-52) 

A(Sjc) a 8(s) =• 0.0 c-3,4.,5* 

From this development we can see the polynomial form 
filter provides the transformation of variables from Zj to 

K a | | for i = i, ,n. After solution of the linear equations 

the filter zeros and poles can be found by factoring the numer- 
ator and denominator. 

3.3 Combined Compensation 

Having considered feedback and cascade compensation 
individually we will now examine the combination of both types 
of compensation. Figure 3-1. Attention is given first to 
tachometer-cascade compensation, Figure 3-8* which indicates 
nonlinearity may be a problem. In this section the system 
gain is assumed known including the filter gain. Naturally, 
the three required equations are from velocity constant and 
root specification. 

The characteristic equation is 

*- (pn + ppi ppiPz+K *■'*)$' 

+ (Hpip AX/tt a 'Ak) 5 + x-t ~ o (3-53 } 

from which we apply the root equations and obtain the non- 
linear matrix equation 
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(3-54) 



A 


8 


C 










o 


8 


0 


A 


D 






— 




K 


8 


O 


F 


| 






r 



where 



A- Kv PI PI 

3 = KKv 

3 = -/C 

5 - (Pl+Pz) cv£<t>, ft) i~ cO,l4>± (?) 

F- F $.,($') (3-55) 

6 ~ -P±Pz w*6(f)-(&+Pz)w*Md-co*h(S) 

H ~ <Fz (t) K 

F - aj A2 dff) +(p± + Pz)w,v3Ff) + d$( f) 
s' - k w H d>, ft) 

r - 



Similar to the bandwidth situation this case lends itself to 
the elimination method of solution such that equation (3-54) 
reduces to a quadratic equation in tachometer gain K-t from 
which the compensator pole and zero may be found. As previ- 
ously the two answers give a choice which should operate under 
the same restrictions previously discussed. 

Finally, we shall consider the combination of accelera- 
tion feedback-cascade compensation. Experience indicates the 



28 



acceleration gain will simplify the equations, but for this 
case we find the compensator zero prevents this complete 
simplification. Using the root and velocity constant specifi- 



cations with the characteristic equation 

S*+ ( PJ+-PZ -f-P-t Hk (P1V2+PPH PP2+MM* it)S 

+ ( PJLPZP i H) S + =■ cD 

we obtain the nonlinear matrix equation 





Vf 


e 


O 


a 


rp " 
lit 




o' 




c 


D 


E 


F 




z 


* 


where 


R 


O 


s 


T 


\p«i_ 


1 ! 


H 
— — 



A = iCvPlPZ 

-K 

c = (Pi t pi) ~t~ oj/ 4 JS) 

~ k4>,($) 

£ -• Kat, c z &(f) 

£ - BO),?4> t (f) 



( 3 - 5 6 ) 



( 3 - 57 ) 



( 3 - 58 ) 



H - - - pi pi -{fit n)(0^4>U f) - lu*4> a ( 

£ ~ ^ < 4 ^-- £/ PZio^ 4 it) i (Pi+P2)0Jn~ L d>z (?) 
s~ K(P^4i,(T) 
r - tenfold 



The system of equations can be reduced by elimination (like 



the tachometer case) to a quadratic in acceleration gain K a 



from which the compensator pole and zero can be found. 



In this section we see the simultaneous equations 



29 



easily lend themselves to solution by elimination but appear 
to be close to linear. This is true if a fourth specification 
is given that doesn*t create non I i neari ti es (like bandwidth). 
The linearization method used throughout this thesis can be 
used for the K-^Z or K a Z term. The situation is somewhat 
easier if a multi section filter is used. There is a nonlinear 
product for each section of filter. Therefore, if the number 
of filter sections is even then specify an additional k/2 
complex roots, linearize and solve. 
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Basic System Block Diagram 
F i gure 3- I 
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Third Order System With Tachometer Feedback 

Figure 3-2 




Third Order System With Acceleration Feedback 

Figure 3-3 
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Third Order System With Tachometer And 
Acceleration feedback 

F i gure 3-4 




Third Order System With Single Cascade Compensation 

Fi gure 3-5 
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Third Order System With Multi section Cascade Compensation 

Factored Form 

Figure 3-6 
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Third Order System With Multisection Cascade Compensation 

Polynomial Form 

Figure 3-7 
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Third Order System With Cascade 
And Tachometer Feedback Compensation 

F i gure 3-8 




Third Order System With Cascade 
And Acceleration Feedback Compensation 

Figure 3-9 
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IV. Examp I es 

Examples in this section were solved on an IBM 360 
computer using subroutine programs furnished with the 
computer. The programs used are specific in nature and 
appear in Appendix II. The following points apply to all 
examp I es : 

a) negative feedback is desired so that negative 
feedback constants are not acceptable since this implies 
positive feedback. 

b) when a cascade pole or zero is found positive 
sign indicates left half plane and negative sign indicates 
ri ght half p lane. 

c) for bandwidth and f eedback-cascade combination 
examples, there are no logic statements in programs to select 
the proper set of parameter values for determining the system 
roots. This part of the program was written with a priori 
knowledge. For roots the sign convention mentioned above does 
not apply, use the standard sign convention. 
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E xamp ! e 4.1 




Using feedback compensation design a system with 
open loop transfer function in Figure 4.1. 

Root S peci fi cation: Zeta = 0.5, Natural Frequency = | o . 0 
No Velocity Constant or Bandwidth Specification 
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INPUT DATA 
**** ****** 

POLES ARE ZERO, 0.500000 01 , 0.(500000 01 

ROOT SPECIFICATIONS ARE ZETA = 0.50000D 00, 

NATURAL FREQUENCY- 0.100000 02 

NO VELOCITY CONSTANT SPECIFIED 

NO BANDWIDTH SPECIFIED 

OUTPUT DATA 
*********** 

FORWARD GAIN 

G.10000E 03 

TACHOMETER GAIN 
C.80000E OC 

COEFFICIENTS OF CHARACTERISTIC EQUATION 
(ASCENDING ORDER! 

C.10000E 03 0.11C00E 03 0.11000F 02 0.10000E 01 

THE SYSTEM ROOTS ARE 

REAL PART 

-0.10000E 01 -0.50000E 01 -0.50000E 01 

IMAGINARY PART 



0.0 -0.86603E 01 0.R6603E 01 
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Examp I e 4,2 




Design an acceleration feedback system with the open 
looo transfer function shown in Figure 4.2. 

Root spec i f i cat i on : Zeta = 0.5, natural frequency = I 1.0 

No velocity constant and bandwidth specification. 

This example indicates an unattainable solution since 
solution gives negative K and Ka and a right half plane 
root. 
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INPUT DATA 
********** 



POLES ARE ZERO, O.IOOOOD 01, 0.50000D 01, 0.60000D 01 

ROOT SPECIFICATIONS ARE ZETA= 0.50000D 00, 

NATURAL FREQUENCY= 0.11000D 02 

NO VELOCITY CCNSTANT SPECIFIED 

NO BANDWIDTH SPECIFIED 

OUTPUT DATA 
*********** 



FORWARD GAIN 
-0. 100 10E OA 

ACCELERATION GAIN 
-0.82645E-01 

COEFFICIENTS OF CHARACTERISTIC EQUATION 
(ASCENDING ORDER) 

-C.10010E OA 0.30000E 02 0.12373E 03 0.12000E 02 

0.10000E 01 

THE SYSTEM ROCTS ARE 



REAL PART 

-0.3A19AE 01 0.2A19AE 01 -0.55000E 01 -0.55000E 01 



IMAGINARY PART 



0.0 



0.0 



-0.95263 E 01 0.95263E 01 
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Examo I e 4.3 




F i gure 4.3 



Redesign of Example 4.2. 

Root Specifications: Zeta = 0.5, Natural Frequency = 12.0 

No Velocity Constant and Bandwidth Specified 
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INPUT DATA 
********** 

POLES ARE ZERO, 0.10000D 01, 0.50000D 01, 0.60000D 01 

ROOT SPECIFICATIONS ARE ZETA= 0.50000D 00, 

NATURAL FREQUENC Y= 0.12000D 02 

NO VELOCITY CONSTANT SPECIFIED 

NO BANDWIDTH SPECIFIED 

OUTPUT DATA 
*********** 

FORWARD GAIN 

0 • 36000E 03 

ACCELERATION GAIN 
0.29306E CC 

COEFFICIENTS OF CHARACTERISTIC EQUATION 
(ASCENDING ORDER ) 

0.36000E 03 0.30000E 02 0.14650E 03 0.12000E 02 

0.10000E 01 

THE SYSTEM ROOTS ARE 
REAL PART 

-0.60000E 01 -0.60000E 01 -0.37068E-07 -0.37068E-07 

IMAGINARY PART 

-0.10392E 02 0.10392E 02 -0.15811E 01 0.15811E 01 
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Examo I e 4.4 




Design a third order system with acceleration and 
tachometer feedback for the open loop system in Figure 4.4. 
Root Specification: Zeta = 0.7, Natural Frequency = 10.0 

Velocity Constant = 6.0 
No Bandwidth Specification 
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INPUT DATA 
********** 

POLES ARE ZERO, 0. IOOOOD 01. 0.20000D 01 

ROOT SPECIFICATIONS ARE ZE T A* 0.70000D 00, 

NATURAL FREQUENC Y= 0. IOOOOD 02 

VELOCITY CONSTANT IS 0.60000D 01 

NO BANDWIDTH SPECIFIED 

OUTPUT DATA 
*********** 

FORWARD GAIN 

0.37500E 0 A 

TACHOMETER GAIN 
0.16613E OC 

ACCELERATION GAIN 

* i ' 

0. 12933E-01 

COMPUTED VELOCITY CONSTANT IS = 0.60000E 01 

COEFFICIENTS OF CHARACTERISTIC EQUATION, 

(ASCENDING ORDER) 

0.37500E 04 0.62500E 03 0.51500E 02 0.10000E 01 

THE SYSTEM ROOTS ARE 
REAL PART 

-0.70000E 01 -0.70000E 01 -0.37500E 02 

IMAGINARY PART 

0.71414E 01 -0.71414E 01 0.0 
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Examp I e 4.5 






— it 




Figure 4.5 



Design a fifth order system with acceleration and 
tachometer feedback for the open !ood system Figure 4.5/ 
Root Specification: Zeta = 0.5, Natural Freguency = 4.0 

Velocity Constant = 0.5 
No Bandwidth Specification 
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INPUT DATA 
********** 



POLES ARE ZERO, 0.20000D 01, 

C.90000D 01, 

ROOT SPECIFICATIONS ARE ZETA= 
NATURAL FREQUENCY 3 0.40000D 01 

VELOCITY CONSTANT IS 0.50000D 

NO BANDWIDTH SPECIFIED 

OUTPUT DATA 
*********** 

FORWARD GAIN 

0.13349E 04 

TACHOMETER GAIN 
0.11909E 01 

ACCELERATION GAIN 
0 • 65497E— 0 1 

COMPUTED VELOCITY CONSTANT IS 3 



0.60000D 01 
0.10000D 02 

0.50000D 00, 
00 



0.50000E 00 
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COEFFICIENTS OF CHARACTERISTIC EQUATION 
(ASCENDING ORDER) 



0.13349E 04 0.26697E 04 0.10354E 04 0.25400E 03 

0.27000E 02 0.1 OOOOE 01 



THE SYSTEM ROOTS ARE 



REAL PART 



-C.63277E OC 
-0.20000E 01 



IMAGINARY PART 



0.0 

0.34641E 01 



-0.11184E 02 



-0.26026E 01 



-0.11184E 02 



0.26026E 01 



-0.20000E 01 



-0 . 3464 IF 01 
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Examp I e 4.6 
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Figure 4.6 



Design a tachometer feedback system for the open 
looo system in Figure 4.6. 

No Root Specification 
Velocity Constant = 0.°l 
Bandwidth Frequency = 4.0 

Systems roots are for number two gain values since system 
g a 5 n one is negative. 
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INPUT OATA 
********$$ 



POLES ARE ZERO, 0.50000E 01, 0.60000E 01 

NO ROOT SPECIFICATION 

VELOCITY CONSTANT I S= 0.91000E 00 

BANDWIDTH FREQUENCY I S= 0.40000E 01 

OUTPUT DATA 
* $ * $ * * * * * * * 

SYSTEM GAIN ONE = -0.31554E 02 

TACHOMETER GAIN ONE= 0.20496E 01 

SYSTEM GAIN TWC= 0.41920E 02 
TACHOMETER GAIN TWO= 0.38325E 00 

COEFFICIENTS OF CHARACTERISTIC EQUATION 
(ASCENDING ORDER) 

0.41920E 02 0.46066E 02 0.11000E 02 

THE SYSTEM ROOTS ARE 

REAL PART 

-0.12317E 01 -0.48841E 01 -0.4R841E 01 

IMAGINARY PART 

0.0 -0.31905E 01 0.31905E 01 
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0.10000E 01 



Example 4.7 




Figure 4.7 



Design an acce I erometer feedback system for the eoen 
loop system Figure 4.7. 

No Root Specification 
Velocity Constant « 5.0 
Bandwidth Frequency ” 3.0 

System roots are for gain values number one since ka 9 
i s negati ve. 



INPUT DATA 
********** 

POLES ARE ZERO, 0.30000E 01, 0.40000E 01 

NO ROOT SPECIFICATION 

VELOCITY CONSTANT IS= 0.50000E 01 

BANOWIDTH FREQUENCY IS= 0.30000E 01 

OUTPUT DATA 
*********** 

SYSTEM GAIN CNE= C.60000E 02 

ACCELERATION GAIN GNE= 0.49388E-01 

SYSTEM GAIN TfcC= 0.60000E 02 

ACCELERATION GAIN TWO= -0.50494E 00 

COEFFICIENTS OF CHARACTERISTIC EQUATION 
(ASCENDING ORDER) 

C.60000E C2 0.12000F 02 0.99633E 01 O.IOOOOE 01 

THE SYSTEM ROOTS ARE 

REAL PART 

-C.29863E OC -0.29863E 00 -0.93660E 01 

IMAGINARY PART 

0.25134E 01 -0.25134E 01 0.0 
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Examp I e 4.8 




Compensate the open loop system in Figure 4.8 using 
a single section cascade compensator. 

Root Specification: Zeta = 0.5, Natural Frequency = 6.0 

Velocity Constant = 4.15 
No Bandwidth Specification 
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INPUT DATA 
********** 



POLES ARE ZERO, 0.50000D 00, 0.700000 01 

ROOT SPECIFICATIONS ARE ZET A= 0.500000 00, 
NATURAL FREQUENCY^ 0.600000 01 

VELOCITY CONSTANT IS 0.41500D 01 

NO BANDWIDTH SPECIFIED 

OUTPUT DATA 
*********** 

SYSTEM GAIN 

0.13756E 04 

COMPENSATOR POLE 

0.37847E 02 

COMPENSATOR ZERO 
0 . 39961 E OC 
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VELOCITY CONSTANT! COMPUTED) 



0 • 4 1 50GE 01 

COEFFICIENTS OF CHARACTERISTIC EQUATION 
(ASCENDING ORDER) 

0.54972E 03 0.15081F 04 0.28735F 03 

0.45347E 02 0.10000E 01 

THE SYSTEM ROCTS ARE 
REAL PART 

-0 . 39200E CC -0.30000E 01 -0.30000E 01 -0.38955E 02 

IMAGINARY PART 

0.0 -0.51962E 01 0.51962E 01 0.0 
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INPUT DATA 
********** 

POLES ARE ZERO, O.IOOOOD 01, 0.15000D 02 

ROOT SPECIFICATIONS ARE ZETA* 0.66500D 00, 
NATURAL FREQUENCY* 0.596000 01 

VELOCITY CONSTANT IS* 0.280000 02 

NO BANDWIDTH SPECIFIED 

ROOT SPECIFICATIONS ARE Z ETA* 0.50600D 00, 

NATURAL FREQUENCY* 0.14660D 02 

OUTPUT DATA 
*********** 

FORWARD GAIN 

0.34637E 04 

FILTER ZEROS ARE 

REAL PART 

-0.51508E 01 -0.91788E 01 

IMAGINARY PART 

0.0 0.0 
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Examp I e 4. 9 
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Figure 4.9 



Design a double section compensator for the ooen 
loop system Figure 4,9, 

Root Specifications: Zeta I = 0.665, Natural Frequency I 

5.96; Zeta 2 = 0.506, Natural Frequency 2 = 14.66 
Velocity Constant = 28.0 
No Bandwidth Specification 



FILTER POLES ARE 



REAL PART 

-0 . 14 107 E 02 -0.14107E 02 

IMAGINARY PART 

-0.13817E 02 0.13817E 02 

COEFFICIENTS CF CHARACTERISTIC EQUATION 
(ASCENDING ORDER ) 

0.I6376E 06 0.55482E 05 0.10125E 05 

0.85632E 03 0.44213E 02 0.10000E 01 

THE SYSTEM ROCTS ARE 
REAL PART 

-0.39634E 01 -0.39634E 01 -0.21451E 02 -0.74180E 01 

-C.74180E 01 

IMAGINARY PART 

-0.44512E 01 0.44512E 01 0.0 -0.12645E 02 

0.12645E 02 
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Examp I e 4.10 




F i gure 4.10 



Design a single section cascade comp'ensator with 
tachometer feedback for open loop system Figure 4.10. 

Root Specification: Zeta = 0.5, Natural Frequency = 7.0 

Velocity Constant = 0.5 
No Bandwidth Specified 

These specifications can not be satisfied. Choice one 
gives negative feedback gain and roots indicated. Choice 
two gives undesirable right half plane compensator. 
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INPUT DATA 
********** 



SYSTEM POLES ARE ZERO,, 0.50000E 01, 0.80000E 01 

ROOT SPECIFICATIONS ARE ZE T A= 0.50000E 00, 

NATURAL FREQUENC Y= 0. 70000E 01 

VELOCITY CONSTANT IS = 0.50000E 00 

NO BANDWIDTH SPECIFIED 

OUTPUT DATA 
*********** 

CHOICE ONE 



TACHOMETER G A I N= 
COMPENSATOR POLE= 
COMPENSATOR ZERO= 


-0. 73364E-01 
0.1 5640E 02 
0.60346E 00 


CHOICE TWO 

TACHOMETER GA I N= 
COMPENSATOR POLE= 
COMPENSATOR ZERO= 


0.60888E 00 
-0.46819E 02 
-0.26924E 01 



COEFFICIENTS OF CHARACTERISTIC EQUATION 
(ASCENDING ORDER) 

0.30173E 03 0.11035E 04 0.20664E 03 0.28640E 02 

0. 10000E 01 

THE SYSTEM ROOTS ARE 
REAL PART 



-0.28840E OC 
IMAGINARY PART 


-0.35000E 01 -0.35000E 01 -0.21352F 02 


o 

• 

o 


-0.60622E 01 0.60622E 01 0.0 
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Examp I e 4.11 




Figure 4,|| 



Redesign Example 4.10 for single section cascade 
compensator with tachometer feedback for open loop system 
F i gure 4.11. 

Root Specifications: Zeta = 0.5, Natural Freguency = 10.0 

Velocity Constant = 0.5 
No Bandwidth Specified 

Data on page 61 shows results for Wn = 10,0 and system 
roots for parameter vaiues choice one. If Wn is changed 
to Wn = |4.0 data on page 62 shows both choices of 
parameter values are satisfactory. 
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INPUT DATA 
********** 



SYSTEM POLES ARE ZERO,, 0.50000E 01, 0.80000E 01 

ROOT SPECIFICATIONS ARE ZETA= 0.50000E 00, 

NATURAL FREQUENCY= 0.10000E 02 

VELOCITY CONSTANT IS = 0.50000E 00 

NO BANDWIDTH SPECIFIED 



OUTPUT DATA 
*********** 



CHOICE ONE 



TACHOMETER GAI N= 0.16103E 00 
COMPENSATOR PCLE= 0.34088E 01 

COMPENSATOR ZEPO= 0.14829E 00 



CHOICE TWO 



TACHOMETER G A I N= 0.12470E 01 

COMPENSATOR POLE= -0.21609E 03 

COMPENSATOR ZERO= -0.22957E 02 



COEFFICIENTS GF CHARACTERISTIC EQUATION 
I ASCEND I NG ORDER) 

0.7A145E 02 0.64829E 03 0.16483F 03 0.16409E 02 

0.10000E 01 

THE SYSTEM ROOTS ARE 
REAL PART 

-0. 11786E 00 -0.62909E 01 -0.50000E 01 -0.50000E 01 



IMAGINARY PART 



0.0 0.0 -0.86603E 01 0.86603E 01 
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INPUT OATA 

********** 



SYSTEM POLES ARE ZERO,, 0.50000E 01, 0.80000E 01 

ROOT SPECIFICATIONS ARE ZETA* 0.50000E 00, 

NATURAL FREQUENCY^ 0.1A000E 02 

VELOCITY CONSTANT IS = 0.50000E 00 

NO BANDWIDTH SPECIFIED 

OUTPUT DATA 
*********** 

CHOICE ONE 



TACHOMETER GAIN= 
COMPENSATOR POLE= 
COMPENSATOR ZERO= 


0.29434E 00 
0.46165E 01 
0.21652E 00 


CHOICE TWO 

TACHOMETER GAI N= 
COMPENSATOR POLE= 
COMPENSATOR ZERO= 


0.16032E 01 
0.43558E 03 
0.87812E 02 



COEFFICIENTS OF CHARACTERISTIC EOUATION 
(ASCENDING ORDER > 

0.10826E 03 0.71652E 03 0.24718E 03 0.17616E 02 

0.10000E 01 

THE SYSTEM ROOTS ARE 
REAL PART 



-0.15980E 00 
IMAGINARY PART 


-0.34564E 01 -0.70001E 01 -0.70001E 01 


o 

• 

o 


0.0 -0.12124E 02 0.12124E 02 
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Examo ! e 4.12 




Figure 4.12 



Design single section cascade compensator with 
acceleration feedback for open loop system Figure 4.12. 
Root Specifications: Zeta = 0.5, Natural Frequency = 5.0 

Velocity Constant = 0.5 
No Bandwidth Specified 

System roots are for choice two parameters 
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INPUT DATA 

********** 



SYSTEM POLES ARE ZERO,, 0.50000E 01, 0.80000E 01 



ROOT SPECIFICATIONS ARE ZETA* 0.50000E 00, 
NATURAL FREQUENCY* 0.50000E 01 



VELOCITY CONSTANT IS * 0.50000E 00 



NO BANDWIDTH SPECIFIED 



OUTPUT DATA 
*********** 



CHOICE ONE 



ACCELERATION GAIN* -0.26303E 00 
COMPENSATOR POLE* 0.20500E-44 

COMPENSATOR ZERO* -0.13047E 02 



CHOICE TWO 



ACCELERATION GAIN* 0.37030E-01 
COMPENSATOR POLE* 0.14807E 02 

COMPENSATOR ZERO* 0.59227E 00 



COEFFICIENTS OF CHARACTERISTIC EQUATION 
(ASCENDING ORDER) 



0.29614E 03 0.10923E 04 0.24345E 03 0.4632 2E 02 

0.10000E 01 

THE SYSTEM ROOTS ARE 
REAL PART 

-0.28868E 00 -0.25000E 01 -0.25000E 01 -0.41033E 02 



IMAGINARY PART 



0.0 -0.43301E 01 0.43301E 01 0.0 
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V 



Genera I Proaram 



The oroblems solved in this thesis were done on an 
individual basis with no attempt made to write a general 
oroqram. In retrosoect it is seen that a aenera I program 
could be written provided the type of specifications to 
be satisfied are selected, the maximum system order Inciudina 
any ■ comoensat i on effects is stated, the desired output 
information and format is known, and the system configuration 
is considered. 

Our prime specification was a complex root specifi- 
cation which is dependent on the form of the characteristic 
eguation coefficients. In all The cases considered the riaht 
hand side of the matrix eguation had the same form. The form 
is a function of system type and order, and has two genera! 
forms 




K O 



(5- I ) 



from eguation (2-3) and 

n 









(5=2) 



from eguation (2-4) where 

Gin, ~ 1*0 




The left hand side coefficient matrix also reduces to a 



convenient form 




( 5 - 4 ) 



where R is an unknown aopearing in the j coefficient 



of the character i st i c equation. If a Darticular coefficient 
of the characteristic equation has no unknowns then the 
appropriate place in the coefficient matrix above is replaced 
by a zero. As noticed for the cascade compensation system 
terms similar to the right hand side are Generated in the 
left hand coefficient matrix multiplying the unknown filter 
Dole or polynomial coefficient except the coefficients are 



C\n t = Sum of- 'nt& Peopocr of- L of- 7H£- f/in=fZ Paces 

c->, - - K-n 

The root equations can be generalized by applying them 
to a characteristic equation with coefficients 



Ak - (Lk^ ~t~ bk.S + C-K ^ 

and then the result can be reduced to a given situation by 



U ^ 

J~ (X k w 



( 5 - 5 ) 



and 



y cuuv 



( 5 - 6 ) 



where 



l< =. v\ •+■ ofeoee. of filtcr. 



Givi = AO 
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reading in the orooer value tor the coefficients 0 ^) and 
The velocity constant also has certain orooerties that 
make it conducive to orogrammi nq . The basic form is 



k v K TT*c- 
ft Pc 



(5-S) 



where n is the number of system ooles and m is the number of 
system zeros. The form above aoDlies for systems with unity 
feedback and acceleration feedback. If tachometer feedback 



is used the velocity constant becomes 



m 



Kv= w 7/ _f‘ 



kX 



(P-O) 



// PC X PPf 7f*i 

The remaining situations are for cascaded systems an r 
the form is now 

K* =. P /T TV /' 



y\ 



(5-10) 



~7TP< 7/ pj 

where k is the number of cascade filter sections. The cascade 
acceleration combination is as above but the tachometer 



cascade combination becomes 

m < 

Kv = jrzs 



jf pc fipj + K fie hj 



(5-11) 



rj T y\ H T.I /' CJ 

In all the above situations if any oole or zero is not in the 
system it is replaced by a one where as removal of tachometer 
feedback imolies = 0. 

The final specification presents a rather te H ious oroblem 
as the general form shows 



67 



/ 






K $ y - -t- ix>) 

<*- / 7 A~4 v t |A ^“ 5 \ 

C &n-i s -t- dn~jS + — + lj + (<i V s%^OL^.i% + 



S=jiO b (5- 12) 



where the a^’s are a function of the unknowns. Since the terms 
are squared it is recommended that bandwidth be set up for the 
most difficult case exoected to be encountered. 

Once the soeci f i cati ons have been oreset in the program 
they can be implemented by reading in aporopriate constants, 
then using logic statements the simultaneous matrix equations 
can be formed. Again as in the bandwidth case the nonlinear 
cases using the elimination method for solution should only be 
set up for the maximum system configuation and order expected 
to be encountered. This is because the computer can process 
numbers but is unable to process unknowns specified by alpha- 
betic characters. Designers may find that for cases of this 
nature it may be beneficial to use an aporopriate iteration 
routine mentioned in the next section. A rough general flow 
chart of major ooerations is shown in Figure 5-1. 
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General Program Flowgraoh 
F i qure 5-1 
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VI. Conclusions and Recommendati ons 



The purpose of this paper was to investigate the feasi- 
bility of using a computer for designing with algebraic methods. 
It has been shown that this can be accomplished, in most cases, 
by the linearizing transformation of variables or manipula- 
tion of system analytic form. This was not true for the 
feedback-cascade compensation cases or when bandwidth speci- 
fication was to be satisfied but these problems could be 
solved by application of a straight elimination process. 

The important advantage of this method is it only restrains 
the system sufficiently to solve for the unknowns. Other 
design attempts tended to specify all root locations in imple- 
menting a solution. 

An area of possible future study is converting other 
system specifications into algebraic eguations and analyzing 

technigues of solution reguired. We have dealt with only 

three specifications in this paper. Can we effectively use 

peak overshoot, steady state error, settling time, resonance 

freguency, or rise time? In analyzing these specifications 

non I i near i t i es may appear that defy solution by elimination, 

therefore, it may require using iteration methods for 

solution like Newton's method or steepest decent, or a 

combination of these. The results obtained by these methods 

are dependent heavily on the initial value guess and to a 

lesser degree the quality of convergence to the final value. 
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An interesting ooint for further investigation would 
be the possibility of deriving another definition for band 
width based on a criteria other than magnitude. If band- 



width can be specified as a comolex number 
linearities are not compounded by squaring 
cation wi I I reduce to a two equations simi 
equations. 



then the non- 
and the specif i 
lar to the root 
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